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In most supersymmetic models, neutralino dark matter particles are predicted to accumulate in 
the Galactic center and annihilate generating, among other products, gamma rays. The EGRET 
experiment has made observations in this region, and is sensitive to gamma rays from 30 MeV to ~30 
GeV. We have used an improved point source analysis including an energy dependent point spread 
function and an unbinned maximum likelihood technique, which has allowed us to significantly lower 
the limits on gamma ray flux from the Galactic center. We find that the present EGRET data can 
limit many supersymmetric models if the density of the Galactic dark matter halo is cuspy or spiked 
toward the Galactic center. We also discuss the ability of GLAST to test these models. 
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I. INTRODUCTION 

Observations by a variety of experiments have revealed 
that a great deal of the mass of our universe is dark and 
cold jlj. Despite this growing body of evidence, we are 
still ignorant of the nature of dark matter. 

One of the most promising dark matter candidates is 
the lightest neutralino in supersymmetric models 0. In 
most models, the lightest supersymmetric particle (LSP) 
is stable by the virtue of R-parity || . Often, this particle 
is a neutralino, x°, the partner of the photon, Z-boson 
and neutral higgs bosons. This candidate is attractive 
due to the fact that it is electrically neutral, not col- 
ored and naturally has the approprite annihilation cross 
section and mass to provide a cosmologically interesting 
relic density. 

Many methods have been proposed to search for ev- 
idence of supersymmetric dark matter. These include 
experiments which hope to measure the recoil of dark 
matter particles elastically scattering off of a detector (di- 
rect searches) , experiments which hope to observe the 
products of dark matter annihilation (indirect searches) 
and, of course, collider experiments Q. Indirect searches 
include searches for neutrinos from the Sun, Earth or 
Galactic center 0, positrons (?J or anti-protons || from 
the Galactic halo and gamma rays from the Galactic cen- 
ter and halo |J. 

Methods of indirect detection which involve the Galac- 
tic center depend strongly on the distribution of dark 
matter in the Galaxy. At this time, there is a great deal of 
debate and speculation over the merits of various Galac- 
tic dark matter halo models. Numerical simulations fa- 
vor models with strong cusps in the central region, such 
as the Navarro, Frenk, White (NFW) and Moore, et. al. 
models |l(J . These models predict increasing dark matter 
density as one approaches the Galactic center, p oc l/r 1 , 
where 7 is 1.0 for the NFW case and 1.5 for the Moore 
case. 



There have been arguments made, based on observa- 
tions, in the favor of flat density core models. These dis- 
tributions, although possible, are probably not capable 
of producing observable signals from dark matter annihi- 
lation, and are not discussed in this letter for this reason. 

Models with strong density spikes at the center of the 
halo have recently receive some attention [jll). In these 
models, cuspy halos generate spikes as a result of adia- 
batic acretion of matter into the central Galactic black 
hole. 

Finally, if halo distributions are clumpy, rather than 
smooth, it would be possible that less dark matter would 
be present in the central region, and the dark matter 
signal diminished. 

In this letter, we will show results for smooth, cuspy, 
halo distributions of the NFW and the Moore profiles. 
These distributions (as well as spikey models) are es- 
pecially interesting to gamma ray experiments, as they 
provide signals from dark matter annihilation which ap- 
pear as point sources. The angular distribution of events 
is proportional to the dark matter density squared in- 
tegrated over the line of sight. A strongly cusped dis- 
tribution produces the vast majority of the annihilation 
signal in an angular region much smaller than the point 
spread function of EGRET, and hence are indistinguish- 
able from a point source. 

II. EGRET POINT SOURCE LOCATION 
ANALYSIS 

EGRET, the Energetic Gamma Ray Experiment Tele- 
scope, launched on the Compton Gamma Ray Observa- 
tory in 1991, is sensitive to gamma rays in the range 
of approximately 30 MeV to 30 GeV. During its oper- 
ation, EGRET's observations included an exposure of 
approximately 2 x 10 9 cm 2 sec in the direction of the 
Galactic center. This paper presents an analysis of only 
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the gamma rays > 1 GeV because the continuum spec- 
trum of gamma rays from neutralino annihilation peaks 
at higher energies and because the point spread function 
of EGRET improves with energy. 

Previous searches for point sources in EGRET data, 
such as the 3EG catalog ]l2| ], used a single mean point 
spread function for each observed gamma-ray above 1 
GeV and spatially binned the data in square bins of sides 
0.5 degrees. However, the EGRET point spread function 
significantly improves with increasing gamma-ray energy 
with 68% of 1 GeV gamma-rays reconstructed within 1.3 
degrees of the true direction as compared to 68% of the 10 
GeV gamma rays within 0.4 degrees. Our analysis uses 
the point spread function as determined by the preflight 
calibration [|l3| for 6 energy bins above 1 GeV, and does 
not degrade the reconstructed gamma-ray direction to 
the nearest 0.5 degree bin. We use a spatially unbinned 
maximum likelihood analysis to determine the best lo- 
calization of a point source. The diffuse Galactic back- 
ground is from the same model used for the produc- 
tion of the 3EG catalog. Before addressing the Galac- 
tic center region, we tested our method on well known 
sources such as the Vela pulsar and the Crab pulsar. 

We found the position with maximum likelihood for the 
EGRET source near the Crab pulsar to be 1=184. 52,b=- 
5.79. The known location of the Crab is 1=184.56, b=- 
5.78 which is within the 95% confidence region as deter- 
mined by of our analysis. The 3EG catalog lists the loca- 
tion of this source as 1=184.53, b=-5.84, with the known 
location well outside of the 95% confidence contour. 

The results for the Vela pulsar are similar. We found 
the maximum likelihood at 1=263.53, b=-2.82 with a 
known location of 1=263.55, b=-2.79 again within the 
95% contour. The 3rd EGRET catalog lists Vela at 
1=263.52, b=-2.86 and the known location is well out- 
side of the 99% confidence contour. 

When our technique is applied to the Galactic center 
region, we find a point source located at 1=0.19, b=-0.08. 
The Galactic center is excluded as the source beyond the 
99.9% confidence level, see figure 1. The 3EG 95% con- 
fidence region includes the Galactic center as shown by 
the circle in figure 1. We find that if this source, modeled 
with a differential power law spectrum with slope deter- 
mined by the maximum likelihood technique to be -2.2, 
is included in the background of the region, the 95% con- 
fidence upper limit on the number of gamma rays from 
a point source at the Galactic center is 10 to 100 (de- 
pending on the spectrum of the source). By contrast, 
the source at 1=0.19, b=-0.08 is a bright EGRET source 
of 370 gamma-rays. The identification of this source is 
unknown |l5| ] , but this new localization agrees well with 
a postulated source of inverse Compton gamma rays from 
the electrons which create the Galactic center radio arc 
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FIG. 1. Unbinned maximum likelihood point source anal- 
ysis of the Galactic center region. 50, 68, 95, 99 and 99.9% 
confidence intervals on the point source position are shown. 
Note that the Galactic center is excluded beyond the 99.9% 
confidence level as the location of the source. The 95% confi- 
dence contour of the 3EG catalog position is shown as a circle 
for comparison. Also shown are all gamma rays above 5 GeV. 



III. CALCULATING THE GAMMA RAY FLUX 
FROM THE GALACTIC CENTER 

We calculated, for a variety of supersymmetric models, 
the number of events EGRET would have been expected 
to have observed, as a function of the halo model. We 
only consider those models which do not violate accel- 
erator limits [plJl^l, including b to S7 18 and invisible 
Z decay width measurements. Furthermore, we require 
that the relic density of the LSP be 0.05 < tt x h 2 < 0.2. 
We calculate the neutralino relic density using the full 
cross section, including all resonances and thresholds, 
and solving the Boltzmann equation numerically [fl9|| . 
Coannihilations with Charginos and Neutralinos are in- 
cluded. We then calculate the LSP annihilation cross 
section, mass and resulting gamma ray spectrum, for a 
given halo model. The results are shown in figures 2 and 
3. 

The general supersymmetric parameter space, even for 
the minimal supersymmetric standard model, consists of 
more than 100 free parameter and, therefore, must be 
simplifed to do any practical calculations. We consid- 
ered a 7-dimensional parameter space consisting of the 
gaugino mass parameter, M2, the non physical mass /z, 
the ratio of higgs vacuum expectation values, tan/3, a uni- 
versal SUSY mass scale, Msusy, the pseudoscalar higgs 
mass ttia, and the couplings A t and Ai,. 

We parametrized the continuum gamma-ray spectrum 
as a function of the LSP mass and calculated the 95% 
exclusion confidence levels which could be placed by the 
EGRET data. The parametrization depends on the neu- 
tralino annihilation branching fractions, but varies lit- 
tle in the majority of models. This exclusion contour is 
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shown in figures 2 and 3 as a solid line. We did not con- 
sider the 77 or 7Z line emission as these are generally 
above the energy range sensitive to EGRET. 

We also considered the ability of the future experiment, 
GLAST, to probe the galactic center for dark matter an- 
nihilations. With larger area and better angular reso- 
lution, GLAST, will be capable of testing many more 
models than EGRET. Furthermore, GLAST, with sen- 
sitivity to energies as high as ~ 1 TeV, can test models 
with heavier LSP neutralinos somewhat more easily than 
EGRET. In figures 2 and 3, the expected sensitivity of 
GLAST, after 3 years of observation, is shown as a dashed 
line. 
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FIG. 2. SUSY model predicted fluxes for a Moore et. al. 
halo profile. Also shown are the 95% confidence upper limit 
of EGRET (solid line) and the expected GLAST sensitivity 
(dashed line). Blue circles represent models with an LSP 
which is more than 95% higgsino, red stars represent models 
with an LSP which is more than 95% gaugino and black x's 
are models with mixed neutralinos. 
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FIG. 3. SUSY model predicted fluxes for a NFW halo 
profile. Also shown are the 95% confidence upper limit of 
EGRET (solid line) and the expected GLAST sensitivity 
(dashed line). Blue circles represent models with an LSP 
which is more than 95% higgsino, red stars represent models 
with an LSP which is more than 95% gaugino and black x's 
are models with mixed neutralinos. 



IV. CONCLUSIONS 

Our analysis of the EGRET data in the Galactic cen- 
ter region indicates an off-center point source, excluded 
beyond 99.9% as the Galactic center. Considering this 
source as background, we found no evidence of a point 
source at the Galactic center and determined the 95% 
confidence upper limits on the flux of gamma rays as a 
function of the WIMP mass. 

We compared these limits to the flux predicted for a 
variety of supersymmetric models and galactic halo mod- 
els. We find that for very cuspy (or spikey) halo models, 
such as the Moore et. al. profile, the majority of vi- 
able supersymmetric models are excluded by our limit. 
We show that the GLAST experiment will have the sen- 
sitivity to further constrain the Galactic halo profile of 
neutralino dark matter. 
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